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Abstract 

During the long course of evolution, nature has learnt how to exploit quantum effects. In fact, 
recent experiments reveal the existence of quantum processes whose coherence extends over un¬ 
expectedly long time and space ranges. In particular, photosynthetic processes in light-harvesting 
complexes display a typical oscillatory dynamics ascribed to quantum coherence. Here, we con¬ 
sider the simple model where a dimer made of two chromophores is strongly coupled with a 
quasi-resonant vibrational mode. We observe the occurrence of wide oscillations of genuine 
quantum correlations, between electronic excitations and the environment, represented by vibra¬ 
tional bosonic modes. Such a quantum dynamics has been unveiled through the calculation of 
the negativity of entanglement and the discord, indicators widely used in quantum information 
for quantifying the resources needed to realize quantum technologies. We also discuss the pos¬ 
sibility of approximating additional weakly-coupled off-resonant vibrational modes, simulating 
the disturbances induced by the rest of the environment, by a single vibrational mode. Within 
this approximation, one can show that the off-resonant bath behaves like a classical source of 
noise. 

Keywords: Quantum effects in biology. Quantum correlations. Open quantum systems 


1. Introduction 

The existence of coherence, caused by the interference of probability amplitude terms, is one 
of the distinctive traits of quantum mechanics. Oscillatory behaviours, ubiquitously observed 
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in quantum systems, are the consequence of such coherent phenomena. The fact that light¬ 
harvesting complexes were experimentally proven to exhibit oscillatory electronic dynamics has 
stimulated a deep debate about the nature of such oscillations, and, consequently, the possible 
role played by quantum mechanics in biologically functional systems HUEISISIEI. Clearly, 
oscillations can also be found in completely classical systems (the appearance of oscillatory 
electronic dynamics within classical models was studied in iiiTi). Thus, it is of the utmost 
importance to give a precise characterization of these oscillations by means of some quantumness 
quantifier, in order to unveil the basic mechanisms adopted by nature, which in future could 
possibly inspire new energy technologies. 

The question of whether coherence in light-harvesting complexes has a quantum or classical 
origin was addressed in a series of recent seminal works HI 13[Mill][12[niHl[IS- In Ref. 
ED, Wilde et al. faced the problem from the macro-realism point of view using the Leggett-Garg 
inequalities to test whether the system dynamics is compatible with classical theories, while in 
Refs. 113 [12 E3 HI US] ED El the presence of entanglement among the electronic degrees of 
freedom was used to assess the genuine quantum character of the whole system. 

As living objects are embedded in their own environment, the dynamical behaviour of elec¬ 
tronic excitations is necessarily influenced by the presence of other degrees of freedom, mainly 
phonons, whose coupling to the system has been identified to be one of the possible causes 
of efficient transport IfTSl [13 |20| |2D [32 l23]l . Stimulated by this observation, in Ref l24l . 
O’Reilly and Olaya-Castro recently investigated non-classical features of the molecular motions 
and phonon environments in a prototype dimer that can be found in light-harvesting antennae of 
cyanobacteria ll25l . cryptophyte algae ll2^l27ll and higher plants Ii28]l23[30l . A characterization 
of quantumness was performed by means of the Mandel Q-parameter and the Glauber-Sudarshan 
quasi-probability P distribution, whose negative regions in phase space are not compatible with 
any classical description of the coupled dynamics. 

Nevertheless, albeit Q-parameter and R-function negativity represent a significant way of 
quantifying the quantumness of a system, they do not catch thoroughly the ultimate quantum¬ 
ness represented by non-classicality of correlations. For instance, the R-function of an Einstein- 
Podolski-Rosen(EPR) pair is positive, in spite of the deep quantum nature of this state ED. 

In this paper, we want to overcome this drawback by directly quantifying the degree of quan¬ 
tumness developed during the coupled exciton-vibration dynamics through the use of entangle¬ 
ment and of quantum discord between system and environment. On the one hand, entangle¬ 
ment is usually viewed as an utterly fragile property at room temperatures, as it can be easily 
destroyed by decoherence E2. Despite these caveats, however, its presence was predicted in 
different biological processes H [12 ED El ED H ED ED- On the other hand, there can exist 
quantum correlations also in the absence of entanglement, as witnessed by many quantum proto¬ 
cols. These quantum correlations are captured by quantum discord, whose definition originates 
from two definitions of the classical mutual information whose quantized versions turn out to 
be nonequivalent 1331 . Quantum discord provides us a criterion of quantumness that is both 
necessary and sufficient, in contrast, for instance, with the use of the Mandel parameter or the 
R-distribution. In the case of entanglement, we will make use of the negativity El, which is a 
sufficient criterion itself, and of a lower bound for the entanglement of formation l35l . The anal¬ 
ysis will be performed by considering the dimer-exciton system both in the presence and in the 
absence of decoherence effects induced by low-energy modes in the phonon environment. It will 
be also interesting to monitor the quantum character of the bath by quantifying the entanglement 
between the bath itself and the system. 
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2. Results 


2.1. Model 

In Ref. m, it was shown that in prototype dimers present in a variety of biological systems, 
efficient vibration-assisted energy transfer in the sub-picosecond timescale and at room temper¬ 
ature can appear. It was also shown that non-classical fluctuations of collective pigment motions 
are dynamically created. Based on these observations, it was suggested that a connection may 
exist between these fluctuations and the high efficiency of the process. The model employed 
consists of an effective dimer coupled to an undamped bosonic mode. Despite its simplicity, it 
captures the essential features of the problem. 

Let us just briefly recall the physical model. A dimer is composed of two chromophores 
whose Hamiltonian reads 


^ + V{o-\(r2 + 0 - 2 ^ 1 ), (1) 

/= 1,2 

where e, is the energy of the excited level of the i-th chromophore and V is the inter-chromophore 
coupling. The operators cr^ (cr~) create (annihilate) an electronic excitation at site i and are 
expressed in terms of Pauli matrices, erf - {cr^±icry)l2. Note that erf erf and (TjCTj are occupation 
number operators with eigenvalues 0, 1, while A = cr|cr-t- crf(T 2 is a constant of motion of Hgi 
with eigenvalues 0, 1, 2. The dimer is strongly coupled to a quantized vibrational mode of 
frequency Wyj^, the phonon Hamiltonian being 

^vib = ^^vib(^l^i + ^^ 2 ), (2) 

where b\ (b:), j - 1,2, are bosonic operators which create (annihilate) one phonon of the vibra- 

.1 

tional mode of the i-th chromophore, so that one can define the corresponding number operators, 
hj = Finally, the electronic excited states interact with their local vibrational environments 
with strength g. For the corresponding interaction Hamiltonian one has 

%-vib cr-(b] + bi). (3) 

!= 1,2 

Combining Eqs. ([T][^ and restricting the dynamics to one electronic excitation, one obtains the 
effective Hamiltonian 

^ex-vib = h- Ver^ - + ^-) + "vib^-^-’ (4) 

where cr^ - erferf^ - erf erf and cr^ = erf erf + erferf and where = {b\ - b^)! V2 is the relative 
displacement phonon mode. A detailed derivation of Eq. Q is given in the appendix. 

2.2. Exciton-vibration correlations 

The quantum character of the bosonic held was analyzed by the authors of Ref. m using 
both the Glauber-Sudarshan P(a)-function and the Mandel factor Q = ((n^) - {h)^)lh - 1. Both 
P(a) < 0 and Q < 0 are used as a sufficient criterion to verify the presence of quantumness in the 
system. Eor instance, a negative value of P(a) implies that the density matrix cannot be expressed 
as a statistical mixture of coherent states and, then, does not admit a classical interpretation. 
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Taking the numerical parameters from the cryptophyte antennae phycoerythrin (PE545) it was 
shown that both Q{t) and the f’(a)-function exhibit partial non-classical behaviour, that is, there 
are regions of time for Q(t) and regions of a for P where quantum fluctuations are detected. 
These regions would vanish in the absence of dipole coupling {V = 0). 

A different approach to nonclassicality has emerged in the literature, which focusses on 
information-theoretic aspects of correlations and has been shown to be deeply inequivalent to 
the quantum phase space and quasi-probability distribution criteria ll^ . Thus, the concepts of 
entanglement and quantum discord, considered as fundamental resources in many quantum in¬ 
formation protocols, offer an independent and conceptually stronger alternative of verifying the 
quantum character of the process under study. 

In order to detect the entanglement between the dimer and the bosonic mode we will resort 
to the Peres-Horodecki criterion, which states that if the density matrix fails to be positive un¬ 
der partial transposition m, then we are necessarily in the presence of entanglement. As a 
quantifier, we use the negativity of entanglement Eff, which amounts to the sum of the nega¬ 
tive eigenvalues of the dynamical density matrix after partial transposition. Indeed, any bipartite 
state has the form g - Cijki\i){j\ ® |k)(/| and the partial transpose map I ®T{g) transforms it 
into g - Yiijki Cijki\i}{j\ ® 10(^1- As any bipartite separable state can be written, by definition of 
separability, as p = PiBa ® ^ ®’ T(p) would map it into pA - p.g^^ ^ )^, which 

is a perfectly acceptable density matrix. Then, all the eigenvalues T, of pA ^re real, positive, 
and obey 2/^/ = 1- If’ on the other hand, some of the eigenvalues of are negative, we can 
conclude that g does not admit a factorized form. Based on these considerations, the quantity we 
are going to calculate is an entanglement witness, even though it is not a proper measure, aside 
from some special cases (a recent example can be found in Ref. iJTll l. States that are entangled 
even though their negativity vanishes are known as bound entangled states, their main charac¬ 
teristic being that it is not possible to obtain pure entangled states from them by means of local 
operations and classical communication. 

The negativity is also useful to calculate a lower bound for the entanglement of formation of 
a bipartite state llTSl . Indeed, entanglement of formation represents one of the most meaning¬ 
ful measures of entanglement, as its regularized version quantifies the minimal cost needed to 
prepare quantum states in terms of EPR pairs. 

Together with entanglement, it is also of interest to monitor the behaviour of quantum dis¬ 
cord £), which is a more general definition of quantumness with respect to entanglement, being 
nonzero even in case of factorized states whose correlations do not admit any classical interpre¬ 
tation ll33l . The definition of discord is given in the appendix. Here we anticipate that, given 
two parties A and B the quantihcation of the correlations between them is the goal of our study, 
it measures the minimum amount of disturbance introduced in the state of party A because of a 
measurement process performed on party B. 

The behaviour of these two quantities is shown in Eig. [^as a function of time. The system 
is prepared at f = 0 in the Gibbs (thermal) state at room temperature T - 270 K of one of the 
vibrational modes, which is initially uncotTelated with the dimer state |2f+): 

= (5) 

Here, //ei|X±) = A±\X^}, where /l± = ± -^(ei - £ 2 )^ + 4-V^/2. We assume that the frequency of 
the mode is much larger than the thermal energy scale Wyj^ » evolution is 

calculated by solving the Liouville-von Neumann equation p = “'Wex-vib’T’]' principle, as 
we deal with an infinite-dimensional system the eigenstates of which are not Gaussian functions, 

4 




t(ps) 


Figure 1: Quantum discord (red, dashed), entanglement negativity (blue, dotdashed) and entanglement of formation 
lower bound (black) as a function of time. The parameters used, taken from Refs. EH ill are the following: ej - S 2 = 
1042 cm“', y = 92 cm“’, = 1111 cm“*, and g = 267.1 cm“*. 


we would need to calculate an infinite number of matrix elements in order to determine the 
exact full dynamics of the state and its correlations. However, taking into account that, in the 
regime tUyib » KbT, at the initial time only a few matrix elements are significantly populated, 
a truncation in the number of excitons is a very good approximation. The truncation scheme 
consists of neglecting all the matrix elements between the threshold h and n + 1. In the model 
under study we can safely take « = 5, as the results do not change (within the machine eiTor) for 
higher thresholds. 

Negativity of entanglement is calculated by applying transposition to the dimer part of the 
density matrix while leaving the bosonic mode unchanged. In the case of quantum discord, the 
dimer represents the part under measurement. In principle, an optimization over a complete 
positive operator valued measure (POVM) with elements {E^} should be performed in order to 
get the optimal measurement. However, for the sake of simplicity, we will limit ourselves to the 
class of orthogonal projectors. In all the cases known in literature, orthogonal projectors give a 
very tight bound and can be safely used without any appreciable qualitative change. 

Working in the Wvib » KbT regime implies that the dynamics is largely dominated by the 
coherent oscillation between |V+, 0) and |X_, 1), which absorbs most of the spectral weight of the 
whole density matrix. We then expect, out of the exact dynamics, a kind of “two-qubit” Rabi- 
like oscillation that unavoidably generates entanglement. As it can be observed, the quantumness 
witnessed by the presence of negativity and quantified by bipartite discord and entanglement of 
formation, is actually present for any f > 0, even in the time windows when the Q-parameter and 
the P-distribution fail to detect it, that is, when the Rabi-like oscillations reach their minimum 
ED. Around those regions, the persistence of quantum effects clearly indicates that, beyond the 
main oscillation, the multi-mode character of the vibration plays an important role. Interestingly, 
Fig. 0 shows that all quantifiers are close to the values corresponding to maximally entangled 
states for most of the time interval analyzed (cf. the value 1/2 for the negativity). We also em¬ 
phasize that, since entanglement is successfully detected by the negativity for any f, the problem 
of determining the possible occurrence of bound entanglement is removed ll^ . 
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2.3. Single-mode description of thermal noise 

So far, we have discussed the interaction of the dimer with a single high-energy vibronic 
mode. In general, low-frequency phononic modes are also present and cannot be neglected, as 
their existence, albeit taken into account at a higher perturbation order, would cause decoherence 
and dissipation of the system under investigation. 

A phonon bath is described by a collection of independent harmonic oscillators that can 
cause incoherent transitions between the system eigenstates. The position of each phonon is 
indeed coupled to the exciton operator cr^, while it is decoupled from the vibrational mode. The 
bath Hamiltonian can be written as 


Hb = ^Ukblbk. (6) 

k 

while the system-bath interaction Hamiltonian takes the form 

^ 8k(o-i ® llvib)(^l + bk). (7) 

k 

In principle, a hierarchical expansion of the interaction could be employed to obtain the 
reduced system dynamics Il38l . The second order (BIoch-Redfield) perturbation theory would 
greatly simplify the calculation. However, due to the strong detuning between the system and 
the bath, non-Markovian effects, which would not be captured by that method, are expected 
to be relevant. Remarkably, a qualitative description of the phenomenon can be obtained by 
drastically simplifying the approach. In the absence of the environment, as illustrated in the 
previous section, the dynamics is deeply influenced by the approximate degeneracy of the levels 
|X+, 0) and |X_, 1). As the initial population of |X+, 0) at room temperature is close to 1, coherent 
oscillations between those two levels are observed which, among other effects, also determine 
the establishment of the quantum correlations described in Fig. [T] As the eigenmodes of the 
bath lie in a region of the energy spectrum that is far apart from the frequencies of the closed 
system, it is natural to ask whether and to what extent the internal structure of the environment 
matters. To this end, we employ a minimal approach replacing the whole environment with a 
single, low-energy bosonic mode Icq. Since, because of the approximation method used, we are 
not in the presence of a true bath, as in any few-body problem, a continuous flow of information 
(which is the cause of oscillations) between the system and the low-frequency mode is expected 
to take place instead of the relaxation behaviour typical of decoherence. Notice that in this way 
the non-Markovian character of the evolution (even if in an approximate form) is kept. 

Using the coupling go and the frequency a>o of ko as free parameters, we explored different 
regimes and found behaviours that are in qualitative agreement with the whole bath case. In Fig. 
I^we plot the population of the excitonic eigenstate \XJ), i.e. Px_ - TinPn,X-(t) as a function 
of time and of go for oiq - lO^^Wvib. In the range of go chosen (centered approximately around 
go - 10 *g), we observe the transition from the coherent regime to the incoherent one, where the 
spectral weight of the |X+, 0) ^ |X_, 1) is reduced. As already pointed out, within our simplified 
model, Fx_ does not reach a true stationary state. However, there is the clear tendency to get 
stabilized around a plateau. As expected, if the frequency of that mode is too close to Wvib, 
where the weak-coupling approximation breaks down, the single mode is not able to capture 
the essential features of the whole bath. This is illustrated in Fig. where we have chosen 
0)0 - 10 *cuvib- In this regime, decoherence is expected to take place before coherent oscillations 
are established El, while we observe high-visibility oscillations. 
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Figure 2: Top panel: Px_ as a function of time and go- The system parameters are the ones given in Fig. ^ and 
cjQ = 10“^ujvib- Lower panel: Px- as a function of time for go = 0 (red) and for go = 10“*g (black). 



Figure 3: Px_ as a function of time and go for tuo = 10 
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Once established the conditions under which modeling the bath as a single mode represents a 
suitable approximation, we use this technique to study the behaviour of entanglement negativity, 
represented in Fig. for the case where decoherence is expected not to completely suppress 
the coherent oscillation (Fig. [T]). For the sake of clarity, here, we will omit the complementary 
discussion about entanglement of formation and discord, as there are no qualitative differences. 
As in the unperturbed case, we observe an initial growth after which the indicator starts going 
down. The important point is that remains positive at any time, showing the resilience of 
entanglement against noise. It may also be interesting to see whether the decrease of dimer- 
vibration entanglement is somewhat compensated by the creation of entanglement between the 
dimer and the bath. Actually, apart from some special values of go and coq and possibly because of 
the roughness of the single-mode approximation performed, this does not happen (see Fig. |5),or, 
at least, entanglement negativity does not reveal it. This means that a global loss of quantumness 
takes place even in the single-mode approximation. Let us stress that, even if we have modelled 
the bath as a genuine quantum system, the results indicate that it actually behaves as a classical 
source of noise. In fact, it induces decoherence in the exciton-vibration part of the model without 
getting quantum correlated itself. Notice that we have considered a very small environment that 
unavoidably presents recurrences regimes and back-flow of information. The absence of bath- 
system quantum correlations is an interesting result, even if not completely unexpected due to 
the fact that the distribution of entanglement among many parties is constrained by the property 
of monogamy and the major quantity of entanglement is already established between the dimer 
and the resonant mode Il39l . 

3. Conclusions 

Understanding the very fundamental mechanism responsible for high-efficiency energy trans¬ 
fer in photosynthesis is expected to lead to both fundamental and practical implications. On the 
one hand, excitation energy distribution is unavoidably influenced by the presence of environ¬ 
mental degrees of freedom, while, on the other hand, the role played by quantum mechanics in 
biological structures has yet to be fully understood. 

We have investigated the presence of quantum correlations in a dimer-exciton system using 
negativity of entanglement and quantum discord. From a qualitative point of view, we have 
found traces of quantumness at any time of interest. This implies that quantum effects are even 
deeper than what indicated by the time behavior of the Mandel parameter or of the Glauber quasi¬ 
probability distribution. Our results suggest that the information conveyed by entanglement and 
discord is richer and could be extended to different models and working regimes. It seems that 
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Figure 5: Negativity of entanglement between the dimer and the low-frequency mode as a function of time and go for 
OJQ = 10 tUvib- 


the phase-space characterization of non-classicality is more sensitive to the main contribution to 
the dynamics, while entanglement and quantum discord also capture the multi-mode structure of 
the bosonic mode. Let us stress that the phase-space parameters, as well as entanglement nega¬ 
tivity (apart from some special case llJTl ') are sufficient criteria for the existence of quantumness, 
while discord is both sufficient and necessary. Then, quantum correlations are established im¬ 
mediately after the interaction takes place. Let us point out that we have observed a correlation 
between transport and non-classical properties. This does not necessarily imply the existence of 
a causal relation, which we would not able to prove, and indeed such kind of answer has not 
been found yet in all the existing literature of the field. Nevertheless, our results point out that 
quantum correlations can play a significant role and may stimulate further studies on the subject, 
eventually addressed to understand whether a functional role exists or not. 

It is also important to set the limits under which non-classicality is robust against thermal 
noise. We proposed a simplified approach to this problem based on the use of a single mode to 
mimic the role of the environment. The computational advantage of this method, compared to 
hierarchical expansions, is evident. As expected, this simplified approach is only meaningful in 
the the weak coupling limit. We tested the method in our system, found the regime where it can be 
applied, and used it to assess the quantum character of the dynamics also in the presence of noise. 
Unexpectedly, despite the high-degree of non-Markovianity and and back-flow of information 
induced by such a small environment, negligible quantum correlations are detected between the 
system and the bath itself. We have found a regime where even a single mode acts as a classical 
source of noise as distribution of entanglement among many parties is subject to monogamy 
restraints. 

Finally, let us stress that the very same model study can also be found in different physical 
contexts. For instance, it could describe electron-phonon interaction in metals. It might also 
be interesting to analyze the dynamical behaviour by considering different system parameters in 
order to establish whether the observed phenomenon is a general characteristic of the model or 
it requires specific experimental condition to be matched. 
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Appendix A 

The model introduced through Eqs. 011 ' is defined in the tensor product Hilbert space i) = 
® i )2 ® ® 5^2, where i )2 are the two-dimensional Hilbert spaces of the chromophores 

and 5^1, 5^2 are the oo-dimensional Fock spaces of the phonons. The analysis simplifies consider¬ 
ably by introducing the eigenstates |A±) of H^\, H^\\X±) - A±\X±), with energy splitting A+ - A- 
- -y^(ei~^-e2)^”+~4f^, and the collective phonon modes = (b[ + b^)/ V2, with b[ (bl) corre¬ 
sponding to the center-of-mass (relative displacement) phonon mode. In view of the properties of 
N, the effective chromophores Hilbert space reduces to the single two-dimensional one-particle 
(or spin j) space with .^i (g> jo 2 '—> Furthermore, since the center-of-mass mode b\ is not 

coupled to the electronic degrees of freedom, only the relative displacement bosonic operator b^_ 
is relevant to the system dynamics, so that the phonon space maps into the single oo-dimensional 
Fock space i.e., '—> i?- Eqs. 03 are combined resorting to the above simplifications 

and hence giving the effective exciton-vibration Hamiltonian introduced in Fq. 0- There, the 
Pauli matrices cr^ and (Tx live in S). 

Appendix B 

Given a bipartite density matrix g and its reduced states = Tr^lp) and g^ = Tr^jp), 
discord can be defined as the measure of how much disturbance is introduced when trying to get 
information about party A when party B is measured ||33. It is defined as 

SAa-.b - - STa-.b, (- 1 ) 

where J(p) = S(gA)-S(gB)-S(g) is the quantum mutual information, obtained from its classical 
counterpart replacing the Shannon entropy with the von Neumann entropy S (g) - -Tr{p log p), 
and where the classical correlations are given by 

Ja-.b = max[5(pA) - 5(A|{£f))], (.2) 

j£fi) 

with the conditional entropy 5(A|{£'^)) = ZjPj‘S(B aie^)’ Pj - ^t^ab^E^q) ^nd where = 
E^glpj is the density matrix after a POVM with elements {E^} has been performed on party B. 
Notice that the definition of discord is not symmetric under the exchange of the two parties. It 
is even possible to find states that are quantum-classical, that is, states that behave as quantum 
objects if one of the parties is observed and as classical objects by observing the other party. 
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